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Abstract 

We compute the momentum broadening of a heavy fundamental charge propagating through 
a M = 4 Yang Mills plasma at large t' Hooft coupling. We do this by expressing the medium 
modification of the probe's density matrix in terms of a Wilson loop averaged over the plasma. 
We then use the AdS/CFT correspondence to evaluate this loop, by identifying the dual semi- 
classical string solution. The calculation introduces the type "1" and type "2" fields of the thermal 
field theory and associates the corresponding sources with the two boundaries of the AdS space 
containing a black hole. The transverse fluctuations of the endpoints of the string determine 
kt = \^XT 3 t: - the mean squared momentum transfer per unit time. (7 is the Lorentz gamma 
factor of the quark.) The result reproduces previous results for the diffusion coefficient of a heavy 
quark. We compare our results with previous AdS/CFT calculations of q. 
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I. INTRODUCTION 



The suppression of high transverse momentum particles in nucleus-nucleus collisions 
known as jet quenching can be used to calibrate the properties of the hot and dense matter 
formed at RHIC [1, 2]. This suppression is the result of the energy loss of partons propagat- 
ing through the medium. For sufficiently high transverse momentum, radiation is expected 
to be the dominant energy loss mechanism. 

There are several approaches to computing the radiative energy loss [3-5]. Here we will 
focus on the multiple soft scattering approximation, known as the BDMPS approach after 
its authors [4, 5]. In this formalism, the medium induced radiation is a function of the 
transport coefficient q = 2kt, which is the mean squared momentum transfer per unit time 
that is imparted to the gluon as it traverses the medium 1 . The value of q was estimated from 
the suppression of high p T particles observed at RHIC, q ~ 10 — 15GeV 2 /fm [6, 7]. This 
value is surprisingly large when compared to the typical momentum scale of the medium 
T ~ 300 MeV. 

The strong jet quenching and the strong elliptic flow [8-11] suggest that the matter is 
strongly interacting. For this reason, it is useful to have a foil to weak coupling calculations 
of transport in high temperature gauge theories. This foil is provided by M = 4 Super Yang 
Mills (SYM) at large t' Hooft coupling, which is tractable once the AdS/CFT conjecture 
is accepted [12-14]. This conjecture states that M = 4 SYM is dual to Type II B string 
theory in an AdS 5 x S 5 background. The correspondence is particularly attractive because 
real time correlators in the strong coupling limit of the gauge theory can be determined by 
solving classical super gravity equations of motion. 

In the field of heavy ion phenomenology, interest in the AdS / CFT correspondence began 
when Policastro, Son and Starinets determined the shear viscosity (77) to entropy (s) ratio 
of the M = 4 SYM plasma [15] 



This ratio says that the mean free path is a fraction of the inverse temperature, £ m fp ~ 
l/(47rT). Since an 77/s ratio of this order is necessary to have hydrodynamics at RHIC, the 
computation in M = 4 SYM was important, because it showed that 77/s can be this small, 
at least in some specific theories. 

Since then there has been a lot of activity in this field. The drag coefficient of a heavy 
quark was computed in Refs. [16, 17] and agrees with the diffusion coefficient found in 
Ref. [18]. The drag also has been computed in different backgrounds [19-21, 23, 24]. The 
properties of bound states of heavy quarks was studied in Refs. [25-30], and the fields 
associated with the jet's passage were studied in Refs. [31, 32]. In addition, the conjecture 
has been used to describe the initial collision [33, 34], and the subsequent hydrodynamic 
evolution [35-37]. 

In Refs. [44, 45], the "jet-quenching parameter" was computed by taking the dipole for- 
mula [38] as a non-perturbative definition of q [4]. The computation proceeds by evaluating 
a light like Wilson loop running along the lightcone. This computation has been extended 
in different ways by several authors [39-42]. 



The more common definition is the mean squared momentum transfer per unit length rather than time. 
The two definitions coincide in the ultra-relativistic limit. 
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In this paper we will define kt as the mean squared transverse momentum transfer to a 
heavy quark propagating through the medium [43]. Along the lines of Ref. [18], we compute 
the medium modifications of the heavy quark density matrix, which in turn is related to 
the momentum broadening. 2 The (1,2) structure of the density matrix is identified with the 
two boundaries of the Kruskal plane of the AdS black hole. The computation proceeds by 
fluctuating an appropriately defined Wilson line. A careful analysis of the fluctuations leads 
to 

Kr = V7AT 3 7r 1 (1.2) 

which diverges in the ultrarelativistic limit. This value is different from that of Refs. [44, 45], 
which is numerically close to the zero velocity limit of this expression. An analysis of the 
approximations underlying the dipole formula used in these works, and the relationship of 
q to the squared momentum transfer, may shed light on the discrepancy, and clarify the 
dynamics of the strongly coupled medium. 



II. TRANSVERSE MOMENTUM BROADENING 

We will for simplicity consider a scalar heavy "quark" coupled to gauge fields. The heavy 
quark is described by the heavy quark effective Lagrangian 

£ = £ YM + (iu- D - M)Q . (2.1) 

Here we follow the conventions, u M = (7, 7V), = (1, v), = 8^ + iA^, and the link from 
the origin to dX is U(dX, 0) = e - idX ^ A » 

Following the formulation of QCD kinetic theory [46, 47], we define the Wigner distribu- 
tion 

f cd (X, r±) = U ca (X, X + r ± /2) Q a (X + vj2)pQ\{X - vj2) U bd (X - r ± /2, X) , (2.2) 

where X denotes the four vector, X = (t, x), r_i_ is a transverse displacement, and X + r± 
denotes (with a small abuse of notation) the space time point, (t, xj_ + r±,z). Then the 
color indices are traced and the Wigner function is averaged with the density matrix of the 
gauge + quark ensemble 

f(X, fl) = (/ CC (X, r ± )> = Tr [p Q\{X - r ± /2)U ab Q b (X + r ± /2)] , (2.3) 

with U ab the straight link U ab (X - r±/2, X + r ± /2) . 

To motivate this definition and subsequent developments, we note that in kinetic theory 
the Fourier transform of this object is identified with the phase space distribution 

f(X, k ± ) = J d 2 r ± e~ ik ^ f(X, r ± ) . (2.4) 



2 The momentum broadening can be expressed as a Wilson loop which is similar to that of Refs. [44, 45], 
but with the Wilson lines approaching the lightcone from below. 
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In this way the transverse current is 

L k ± /(X,k ± ) = -zA /(X ,r ± )=Tr 



p {Q'[--b ± q +\--n ± Q) q 



(2tt) 2 ' ^ dr ± - 

'(2.5) 

With this introduction, we will identify the mean transverse momentum squared with 

(pi(*)> = jd^J |^ k i/( X ' k J' ( 2 - 6 ) 

= J d*x-V 2 r J(X,r ± )\ r±=0 . (2.7) 

If the heavy particle starts with a narrow transverse momentum distribution, then after 
a long time T we expect that the average squared transverse momentum is 

<pi(T)) = 2 KT T, (2.8) 

where the factor of two accounts for the two transverse directions, and kt is the momentum 
diffusion coefficient. 

We next study the time evolution of the Wigner function. A complete set of states with 
one heavy quark may be written with the short hand notation 

J2 I Q\l)\A 1 ){A 1 \Q{\)= [ rf 3 x 1 gl 1 (x 1 )|^)(A M |g Ql (x 1 ), (2.9) 

where A\ labels the Eigenstates of the gauge fixed operator A M (x) (see Ref. [48] for a convinc- 
ing discussion of constraints). Then, working in the Schrodinger picture, we first evaluate 
the density matrix 

Tr[p(t)] = Tr [ e - m ^ p(t ) e +iH ^ ] , (2.10) 
were p(t Q ) is the density matrix at time t a . Inserting complete set of states we obtain 

Tr[p(t)] = \ (AolQme-^-^Q^DlA,) 

a oAi a 2 J te£ 

x WQiDpiQQ^^lAz) 

x (A 2 \Q{2)e +iH ^Q\0)\A ). (2.11) 
Putting the first term at the end, we rewrite this as the path integral 

[ j [DA^] [DQDQ ] ] p° aia2 [x l7 x 2 , A x , A 2 ] e^c **oCyu«**»d-mq ^ t _ ie) Q ^ 

(2.12) 

with 

p° l02 [xi,x 2j ^i,A 2 ] = (AtlQa^MQQK^lAi) . (2.13) 

Here the path integral is performed along the closed time path, starting from t Q , running up 
to time t, and returning back to t Q — ie. 

The energetic quark may be integrated out of this expression for the density matrix. The 
contour Green Function of the heavy quark field in a fixed gauge background is written 

zC(2,l) = (T c Q a2 (x 2 t 2c )Qt i(xi t lc )) , (2.14) 



X Q — (t , %o) tc 

* 9 ^ X = (t, x + v At) 

Xf = (t — ie, x — ie) t' c 

FIG. 1: The contour world line of the heavy quark with velocity v in the z direction. The world 
line of the quark passing through X Q may be parametrized as Xx o {t c ) = Xtf + v^ifo — t oC ), with 
= (l,v). The Wilson line W c [0] follows the world line of the quark. The circles at t c and t' c 
indicate the insertions of the field strengths F yfJ, v^(t c ) and F yi, v u (t' c ) into the Wilson line as in 
Eq. (2.25) . 

where a, denote color indices, tjc is the contour time, and T c denotes the contour ordered 
product. This Green function satisfies 

[iu-D- M) iG{2, 1) = z<W 3 (x 2 - xi)$ (t 20 - ti c ) , (2.15) 

which has solution 

<G(2,1) = e + iMu < x *- x ^ [ ^6(t c -t lc )6UX 2 -X Xl (t c )) 

Jc 7 

Pexp (-Z [ C dt f c v^{X Xl {Q) 

V J t 1C 



X 



(2.16) 



a2ii 



where X^^o) = + v^ L {tc ~ tic) is the heavy quark world line which passes through X\. 
Then integrating the quark fields we have 

Tr[p(t)]= I /"[DA M ]e^o^^Mdet(m-D-M)p° ia Jx o ,x o ,A,A 2 ]W c [0] aia2 , (2.17) 

J x J 



where 



W c [0] = T c exp -i / dt c v^Afj, (X Xo (t c )) . (2.18) 



C 



The path of the Wilson line runs along the world line of the heavy quark and returns on the 
time reversed path as shown in Fig. 1 . 

We make three comments: (i) Up terms suppressed by powers of T/M the determinant 
may be dropped in this expression. This determinant is responsible for effective vertices 
which appear in the light quark Lagrangian when the heavy quark is integrated out [49]. 
(ii) Due to the overall translational invariance of the problem, the average is independent 
of x , and yields an overall factor of volume, (iii) We will denote the contour path integral 
over the gauge fields in Eq. (2.17) as (. . .) A , with the initial and final values of the gauge 
fields given by A\ and A 2 respectively. A\ and A 2 are integrated over with weights given by 
the initial density matrix. Thus Eq. (2.17) becomes 

Tr\p{t)\ = V (trp°[x ,x , A 1 ,A 2 ] W C [0]) A , (2.19) 

where the trace is over the color indices. 
Next we consider the evolution of f(X,r±) 

f(X, rj = Tr [ e -^-M p(t ) e +iH ^ Q\{X - r ± /2)U ab Q b (X + r ± /2)] . (2.20) 
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(t , x ± + t ± /2, z ) (t, x ± + r ± /2, z + uAt) 



(t - ie, x ± - r ± /2, z - ie) (t, x ± - r ± /2, z + uAt) 

FIG. 2: Graphical representation of Eq. (2.23). The Wilson line indicated by the black line is 
denoted FF c [r_i_/2, — rj_/2]. This Wilson line is traced with the initial density matrix, p° aia2 - 

Inserting complete set of states as before we obtain 

f(X, r± ) = J2 I < A °I e ~ m ~ to) ^(i) \ A i) 

A A A A J 0.1.2 



AoAiA 2 A 3 



x (Ax\Q(X)p(t )Q%)\A 2 ) 

x (A 2 \ Q(2)e +lH ^ Q[(X - vj2)U ab \A 3 ) 

x (A 3 \Q b (X + r ± /2)Q\0)\A ). (2.21) 



Putting the first term at the end, we rewrite this as the path integral 

J [DA^] [DQDQ ] ] p° ai Jx 1 ,x 2 ,A 1 ,A 2 ] e^c ^ c c ym+ q^u.d-m)q 



XlX 2 



x Q a2 (x 2 , t - ie) Ql(X - r ± /2)U ab Q b (X + r ± /2) ( Xl , Q . (2.22) 

Performing the quark integration, we have 

f(X,r ± ) = (trp°[x + r J ./2,x -r ± /2,A 1 ,A 2 ] W c [r±/2,-r x /2]> A , (2.23) 

where the Wilson line W c [r±/2, —r±/2] is shown by the black lines in Fig. 2. We may 
expand this for small r_i_ by inserting unity U (X + r±/2, X )U(X , X Q + r±/2) and U (X Q — 
rj_/2, X )U(X , X — rj_/2) at the beginning and ends of the contour and using the definition 
of f ab . We have 

j.2 I Q2 \ 

f(X, rj = f(X, 0) + ^ / ab (X D , r ± )W c [0] ab J + momentum broadening , (2.24) 

with the momentum broadening given by 

1(y) 2 S c S c dtcdt ' c (^P^^A^A^TciF^itc^Fy^QvvWclQ]]^ . (2.25) 

We again make three comments: (i) The first term is in Eq. (2.24) is characteristic of the 
initial momentum distribution of the quark which may be supposed small; the momentum 
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broadening is described by the second term, (ii) We have tacitly assumed that Tr[p(£)] was 
unity, we will divide by Eq. (2.19) to set the normalization, (iii) The T c ordered product of 
field strengths means that we insert field operators into the Wilson line W^cfO] as shown in 
Fig. 1. As in Ref. [18], we define the contour Wilson line W"c[fo/], with deformations in the 
y direction at various points along the contour, 8y(t c ). Then, the contour ordered product 
of field strengths can be written as the variation of W c [Sy] at times t c and t' c respectively. 
The shape of the Wilson line is the source for the contour ordered product of fields. 
Thus from Eq. (2.7), Eq. (2.8), and Eq. (2.25), the momentum broadening is 



- T - \ <^kw A I L dt ° dt '° ( tr <»°- x °- A - smm>) A • (2 ' 26) 

where T is the total real time of the process. 

Now we brake up the contour integration into type "1" and type "2" fields; we write 
5yx(t) and Sy 2 (t') for variations on the one and two branches and use t,t' for the real part 
of t c ,t' c . Then the contour correlations are 

W)= , ' , UP^M) , (2.27, 



(tvp°Wc[0,0]) A \ r 5y 1 (t)5y 1 (f) 

' L-^M), „ 28) 



(tY P °Wc[0,0]) A \ r S yi (t)5y 2 (f) 



{tT P °W c [0,0])A \ S yi (t)5y 2 (t') 
n u 4i\ 1 /+ o ^W C [Sy2,Sy 2 } \ 

iG 2 i(t,t) = „ T rn ( trp - ) . (2.30) 



(tY P °W c [0,0]) A \ r 5y 1 (t f )5y 2 (t) 
Finally, we may use approximate translational invariance to write 



k t = lim - / dte +iuJt (iGnit, 0) + iG 22 {t, 0) + iG 12 (t, 0) + iG 21 (t, 0)) . (2.31) 
ui— *o 4 J 



III. THE ADS/CFT CORRESPONDENCE 

In the previous section we have identified the source for the contour ordered electric 
fields as the variation of the trajectory of a Wilson line running along the contour shown in 
Fig. 1. The strategy to extract the transverse momentum diffusion coefficient kt parallels 
Ref. [18] closely. We first construct a semi-classical string which is the gravity dual of the 
unperturbed Wilson line running along the Schwinger-Keyldish contour. Then we vary the 
endpoint of the string in order to obtain the appropriate contour correlation functions of 
electric fields or Kt- 

A. Preliminaries 

First we recall previous work by us and others to establish notation [16-18, 50] 
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The metric Gmn which corresponds to TV = 4 Super Yang Mills at finite temperature 
is the AdS space with a black hole 

where /(r) = 1 — (y) 4 , i? is the AdS radius, and r is related to the Hawking tem- 
perature iiTR 2 = r . We define the scaled units, nTt = i, nTx = x and f = r/r Q , and 
define z = 1/f , so that the metric reads 

and we distinguish the function / by its argument, f(z) = 1 — z 4 . Further, we will 
sometimes change variables to u = z 2 . In what follows we will drop the "bar". 

The gravity dual of a Wilson line moving along the unperturbed trajectory £3 = vt is a 
semi-classical string stretching in AdSs x S5 whose endpoint follows the curve £3 = vt. 
The dynamics of the semi-classical string is described by the Nambu-Goto action 

1 

Vng - 2 

where 



Sng = / drda\/-deth ab , (3.3) 



h ab = G MN d a X M d b X N . (3.4) 
The solution to the classical equations of motion is described by the map 
(r, a) t— > (t = r, xj_ = 0, x 3 (t, a), z = a, Q 5 = Const) , 

with 

v 

Xs(t, z) — vt-\ — [arctan(z) — arctanh(z)] (3-5) 
B. The Semi-Classical String in the Kruskal plane 

The semi-classical string solution given in Eq. (3.5) is dual to a quark moving with 
velocity v. In the gauge theory this quark is constructed by turning on a U(l) electric field 
and allowing the distribution to come to stationary state [16]. Far in the past the heavy 
quark is nearly at rest, and the world sheet covers the full Kruskal plane [18] 3 . This is 



Our Kruskal conventions are the following. We first define 



= / -]h = ^tan- 1 (z) + itanh \z) (3.6) 

Then we define the coordinates {y+, vJ) 

v + = t + z*(z) v- =t — z*{z) (3-7) 

that for instance f_ = Const is the geodesic of an in-falling lightlike particle. Finally the Kruskal 
coordinates are 

U = -e- 2u + V = e +2v - (3.8) 
We also define V = y e 2tan ~ 1( » and X' = -UVe 2tan '^ z K 
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v=o 



11=0 



V<0 
U>0 




FIG. 3: Kruskal diagram for the AdS black hole. The coordinates (t, r) span the right quadrant. 
The thick hyperbolas on the sides of the two quadrants are the boundaries at r = oo. The 
boundary in the right and left quadrants correspond to the "1" and "2" axes of the thermal field 
theory respectively. The static quark corresponds to a string which spans the full Kruskal plane. 
At finite velocity, the asymptotic solution Eq. (3.5) is discontinuous at t = — oo or V = 0. 



shown in Fig. 3 which also illustrates the Kruskal coordinates. As the electric field is turned 
on in the gauge theory, the quark accelerates and slowly approaches the stationary velocity 
distribution. In the gravity dual this corresponds to slowly turning on a U(l) electric field 
in the boundary brane and waiting for the string to reach the asymptotic form given by 
Eq. (3.5). We re- write this solution in Kruskal coordinates as 



Examining this solution we see a logarithmic divergence for V — > 0, i.e. in the distant 
past, t — > — oo. This discontinuity in the distant past reflects the fact that the asymptotic 
solution Eq. (3.9) is not a solution as the electric field was turned on slowly. At t — > —oo 
it is reasonable expect that this solution slowly deforms into the static solution described 
above and covers the full Kruskal plane. 

Therefore, we will demand analyticity across the V = line while extending V through 
the lower half plane to negative real argument 4 . Through this process x% becomes imaginary 



This imaginary value for the coordinate x% is appropriate. The left quadrant of the Kruskal 
plane is associated with the "2" branch of the field theory with a = i/3/2, or iir/2 for our 



4 This choice of extending V through the lower half plane is compatible with the construction of Son and 
Herzog for the real time path integral which we adopt later. 



£3 = — \og(V) + v arctan(z) for V > . 



(3.9) 




(3.10) 
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rescaled units. Since the operators on the real axis are evaluated at the space time point 
(t,xs) = (t,vt), when time is extended in the imaginary time direction we expect that the 
operators should be evaluated at 

(t,vt)^(t-i/3/2,vt-ivl3/2). (3.11) 

Thus the analytic extension of the string solution into the lower half V plane maps to the 
appropriate X3 coordinate. 

In what follows we will take this string solution as the gravity dual of the heavy quark 
(or Wilson line) propagating along the Schwinger-Keyldish contour shown in Fig. 1 with the 
second axis displaced by if}/ 2. We will make the correspondence between the "1" and "2" 
axes of the contour and the endpoints of the string on the right and left boundaries of AdS. 
When studying fluctuations of this semi-classical string solution, we will bear in mind this 
discussion of the past infinity and demand analyticity across the V = line. 



C. Fluctuations and World Sheet Black Hole 

Having constructed the string solution corresponding to a heavy quark propagating along 
the Schwinger Keldysh contour, we will proceed to fluctuate the shape of the string solu- 
tions in the transverse directions and solve for the fluctuations. Once the classical solution 
is determined in the Kruskal plane we follow the general philosophy of the AdS / CFT corre- 
spondence and equate the classical action of the source to the generating functional 

is NG [SyuSy 2 ] = — L— - (tr p°W[S yi ,6y 2 }) A . (3.12) 



e iS NG [o,o] (trp°W[0 

However in order to find the classical solution, boundary conditions are needed. The appro- 
priate boundary conditions for Gn and G22 etc. are not obvious. 

The appropriate boundary conditions are greatly clarified by a change of coordinates 
which transforms the original induced metric to a diagonal metric which turns out to be 
equivalent to a world sheet black hole. Using the original coordinates (r, a) = (t, z) the 
induced metric is 

hu = -5(^ 4 ) ' (3 ' 13) 
R 2 1 / z A \ 

h u = -~v 2 *>. (3.15) 
z J 

This non-diagonal world sheet metric can be diagonalized by the following change of coor- 
dinates 

t = — — ( t + - arctan(z) — -■y/j arct an (y/^z) — -arctanh(z) + - v /7artanh( A /7z) |3,16) 
y7 \ 2 2 2 2 j 

z = y/jZ. (3.17) 
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and the metric takes the form 



ha = ~*m, (3-18) 

z z 

= ^7^. (3-19) 
z 2 f(z) 

htt = 0, (3.20) 

with f(z) = 1 — z A . Thus the world sheet metric is in fact a black hole with the horizon at 
z 2 = I/7 or z = 1 5 . As usual, the horizon is a coordinate singularity and is a consequence 
the coordinate transformation in Eq. (3.16). In this coordinate system, the induced metric 
for the moving string with "hat" variables equals the induced metric of the static string with 
"un-hatted" variables. 

Next we show how different regions of the world-sheet black hole map into space time. 
First we define the analogs of Kruskal variables for hatted variables z)_, i> + , V, U, z*, . . . . (For 
instance z* = (tan - 1(1) +tanh~ l(z))/2 ). Taking V positive, we have the following relation 
between the world-sheet and space-time Kruskal variables 6 

z = Viz, (3.21) 
V' = (V'f 1 ^ , (3.22) 

X> = — X' = — . (3.23) 

1+z l+z v ; 

From the form of this map we see the upper half of the world sheet Kruskal plane (U, V) gets 
mapped into the upper half Kruskal plane of space-time V > 0. The structure is illustrated 
in Fig. 4. As V becomes negative V becomes imaginary. This is because the solution given 
in Eq. (3.5) is not valid at past infinity. In the lower half (V < 0) of the Kruskal plane we 
can find a separate string solution which maps the lower half (V < 0) of the world sheet 
Kruskal plane. At past infinity, before the electric field was turned on, the quark is at 
rest, and the two string solution are joined. We will therefore connect the solutions in the 
upper half and lower half Kruskal planes by demanding analyticity across the V = line, 
extending V through the lower half complex V plane. This extension through the lower 
complex V plane is consistent with the Herzog-Son construction [52] and physically says 
that only positive energy solutions emerge from past infinity. Again by analogy with the 
Herzog-Son construction, we will extend through the upper complex U plane. Physically 
this says that only negative energy solutions emerge from the future event horizon of the 
world sheet black hole. 

Having clarified the analytic structure of the world sheet black hole and determined 
the appropriate boundary conditions, we next analyze small fluctuations. The equation of 
motion for small transverse fluctuations can be found easily in "hatted" coordinate system. 
Introducing 

y = Viy, (3-24) 



5 This fact has been already noted in [51] 

6 As noted before, we define V = Ve 2tan and X' = —UVe 2tan ^ and analogous relations for V' and 
X' 
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FIG. 4: (a) The world sheet black hole, (b) The Schwarzschild black hole in space time. The 
future event horizon of the world sheet and its corresponding image z = l/y 7 ^ are indicated by the 
solid red line. Similarly the solid blue line maps to the future event horizon of the space time black 
hole. The shaded region of the world sheet black hole maps into the shaded region of space-time 
region of the Schwarzschild geometry. Similarly there is a separate string solution in the V < 
half of the (U, V) plane which maps to the V < half of the (U, V) plane. The two solutions are 
joined at past infinity (V = 0) by demanding analyticity in the lower complex V plane, i.e. that 
only positive energy solutions emerge from past infinity. 



the action for small fluctuation is the same as for the static string [18 

R 2 f didu 



S 



NG 



2na' .1 2m 3 / 2 

We define, 




-4f(u)u(y' 



(3.25) 



y(t,u) = J e- M y(u)Y(u,u)^-. (3.26) 

where we choose to normalize Y&(u = 0) = 1, i.e. y(uj) is the value of the fluctuation at the 
boundary. 

The Euler-Lagrange equation for the small string fluctuations are 

s^-Q+maA + ^i-^o. (3.27) 

This equation is solved by 

Y(u, u) = (1 - up" /4 F(u, u) (3.28) 

where F(u) is a regular function of u > . (1 - u)-™/ 4 is in -falling in the world-sheet 
horizon u — 1. The complex conjugate of this expression is also a solution of the differential 
equation and is outgoing at the horizon. 
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Let us now extend these solutions into the full Kruskal plane. It is useful to express 
the solution in terms of the world sheet Kruskal coordinates. Close to z = 1 the in-falling 
solution and out-going solutions behave as 

e~ lU Y(u;,u) ~ e -^/ 2 MV) i n _f a lli n g (3.29) 
e -v*Y*(u>,u) ~ e^VM-u) ou t-going (3.30) 

As mentioned before, the coordinate V does not cross zero in this z = \fu = 1/ yPy. However, 
U does and we should be careful when defining the branch cut for the logarithm. In analogy 
with the Herzog-Son prescription, we will impose that the solution should be analytic in the 
upper half of the complex U plane. Since U is negative in this quadrant, the prescription 
means that 

e -^Y*(u,u) ~ e *tf/2 e <a/2in(ff) _ (3. 31 ) 

Thus, when the fluctuation crosses the world-sheet event horizon in the right quadrant it 
picks a factor exp{7ru>/2}. Due to the different relation between U and the local (t,z) 
coordinates, in the left quadrant the solution behaves as 

e~ i&i Y*(uj,u) ~ e ^ /21n(f>) . (3.32) 

Thus, after passing the world sheet horizon, keeping the same analyticity properties in U, 
the solution picks a factor exp{— tiuj/2} in the L-quadrant. 

We now address how the matching of left and right quadrant is performed. Using the 
change of coordinates Eq. (3.16) we obtain (for u < 1) 

e~ iU Y{u, u) = e~ iw \\ - u)- iuj/4 Fi(cu, u) , (3.33) 

e -i^Y*(u),u) = e- iut (--u)^ /2 (l-uy iuj/i F {u,u), (3.34) 

7 

where we have defined to = and Fi(u, u), F (u, u) are two regular solutions in u. The 

pole at u = 1/ Ai /7 in the outgoing solution is a consequence of crossing the U = line. Note 
that this pole does not appear in the in-falling solution, since we do not cross the V = 
line. From this point of view, the previous prescription for the analytic properties of the 
solution in the U coordinate translates into the prescription for going around the pole in 
u = l/-k/7- Taking this prescription into account, in the R quadrant the two solutions close 
to the horizon behave as 

e-^Y{u,u) ~ e~ luJ /21n(y) , (3.35) 

Note that from the point of view of the AdS black-hole both solutions are in-falling. The 
exponential factor in the outgoing solution is a consequence of the prescription to cross the 
pole at u — 1. In the same way, the fluctuations in the L-quadrant behave as 

e-^Y(u),u) ~ e -**/*H-V) } (3.37) 
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With these four expression we find four different solutions defined in both quadrants of 
the (AdS) Kruskal plane 



j e~ ia>t Y(Cj, u) in R „ JO in R 

1 m L e Y(uj,u) m L 



j e in R . JO in R 



Following the Herzog-Son prescription [52], we look for linear combinations of these ex- 
pressions that, close to the horizon, are analytic in the lower half of the complex V plane. 
With this criterium only two linear combinations can be found: 

Vo = Vr,o + a yL, , (3-41) 
Vi = VR,i + ctiU Li i . (3.42) 

From the close to horizon behaviors of the solutions Eq. (3.35), Eq. (3.36), Eq. (3.37), 
Eq. (3.38), the analyticity properties demand 

a Q = e +™ e -™/2 j (3 43) 

on = e~™ /2 . (3.44) 

These two solutions are used as a basis for the linearized string fluctuations defined over 
the full (AdS) Kruskal plane 

(a(u)y (u) + b{u))yi{u))) . (3.45) 

This prescription recovers the results of Herzog and Son [52] when 7 = 1. To see this one 
has to realize that when 7 = 1, since these solutions are only defined for u < 1 we do not 
cross the pole and, thus, the exponential prefactors in y{L,R} do not appear. 

The coefficients a(uj),b(u) can be determined by the boundary values of the solutions. 
Thus, if we have 7 

/duj 
—e-^%{u), (3.46) 

y(iu = 0)\ L = I ^Le-**jk{&) , (3.47) 



we obtain 



a(uj 



1 -(-^)+^) e ™/ 2 ) , (3.48) 
b ^ = i^y^) ~ ^ /2 y^)) • (3-49) 



1 



7 Note that at u = 0, t = t/^fy. Thus y(t,u= 0) = / f^e _i<at "#(&), i. e. y(u>) = y(u) 
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We can now compute the boundary action in terms of the string solutions. In (t, u) 
coordinates 



S 



B 



R 2 
lixa 



duo 1 



IJR 



27TM 1 / 



T -y(-u,u)duy(uj,u) 



(3.50) 



Using Eq. 3.48 this action can be expressed as 



Sb = 2^ 



2^ 



(3.51) 



yi (-u) yi (uj) I (h + 1)^Y*(-lj, u)duY(u, u) - h-^Y(-u, u)duY*(u, u) j + 



yi (-u)y 2 (u)e^ 2 n 



u 



\pY\-u,u)d*Y{u,u) + ^Y(-u,u)d & Y*(u,u) j + 



y 2 (-io) yi (u)e-^ 2 (l + n) (-^Y*{-u,u)daY{u,u) + ^Y{-(b,u)d^Y*{Cb,u) j + 



y 2 (-uj)y 2 (uj) (^ : ^Y*(-u,u)duY(ij,u) - (n + l)^^Y(-u,u)d^Y*(u,u) 

where h = l/(exp{7u2>} — 1) and we have used the fact that yx^Cj) = Vi,2(u). From this 
expression we can read off immediately the different correlators by taking derivatives with 
respect to the boundary values. In particular 



lim iG 12 (u) = lim iG 21 (cu) = lim - (iG n (uj) + iG 22 (u)) , 

w >0 U! >0 U) >0 2 

R 2 2 



uj — >o na Tiui 



lim r\/7 Im { tt ^Y*(-u,u)duY(u;,u 



it 



(3.52) 



Since the action for the fluctuation Eq. (3.25) is formally the same as that of the static 
string, we obtain the same result for the transverse momentum transferred as in the static 
case (up to an overall factor of ^/t)- After restoring physical units and using the result of 
[18] and Eq. (2.31) 



Kt = a/tAT 3 7t . 



(3.53) 



IV. CONCLUSIONS 



In this paper we have studied the medium modifications of a heavy quark that propagates 
in a strongly coupled plasma at finite velocity v. Starting from the density matrix of the 
quark, we have expressed the transverse momentum broadening of the probe as a Wilson 
loop running along the x = vt line, with a transverse separation, Ax. This Wilson loop is 
similar to that considered in Refs. [44, 45] but approaches the lightcone from below. The 
second derivative of this Wilson loop with respect to the transverse separation Ax yields kt 



This expression looks slightly different that Eq. (28) of [52] because in that work fk{u) is defined as 
outgoing, while here Y(uj, u) is infalling. This notation is chosen to connect with our previous work [18]. 
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- the mean squared transverse momentum transfer per unit time acquired by the probe in 
its propagation through the medium. We have stressed that the appropriate Wilson loop 
has a time ordered line and an anti-time ordered line. These two pieces can be understood 
as the quark's amplitude and complex conjugate amplitude respectively. Due to this specific 
time ordering, we have introduced the type "1" and type "2" fields of Schwinger-Keyldish 
formalism [53, 54]. The sources for these type "1" and type "2" fields correspond o the 
boundary values of SUGRA fields in right and left quadrants of the Kruskal plane. 

To construct the gravity dual of this ("1","2") Wilson loop we have extended the moving 
quark string solution found in Refs. [16, 17] to the left Kruskal quadrant. The string world 
sheet is disconnected along the V = line because the V = line represents the distant 
past. In the distant past the moving quark string is not a solution to the equations of 
motion because the electric field used to accelerate the quark is being slowly turned on. At 
past infinity the two disconnected solutions are joined by demanding analyticity in the lower 
half plane V as required by the correspondence between AdS black holes and the real time 
thermal field theory. 

Subsequently we have computed by fluctuating the transverse position of the string 
endpoint and solving for the standing waves. To determine the correct boundary conditions, 
we first noticed that the induced metric of the string worldsheet is that of a black hole. 
The event horizon of this world sheet black hole maps to the line r = r ^Py, with r the 
event horizon of the AdS black hole. Then, borrowing from the work of Herzog and Son 
[52] , we have analytically continued the fluctuations through both the world-sheet and AdS 
horizons. This analytic extension amounts to specifying the boundary conditions on the 
physically correct solutions in the full Kruskal plane. From the solution in the full Kruskal 
plane we obtain 

K T = y^XT 3 7T, (4.1) 

which diverges in the v — > 1 limit. 

Finally we wish to compare k,t (the mean squared transverse momentum transfer to a 
heavy quark per unit time) to q. Up to a factor of two, the definition of k,t given here is 
often ascribed to q, at least if the qualifier "heavy" is removed from heavy quark. 9 However, 
the value 2kt differs from the value of the jet quenching parameter found by Liu, Rajagopal 
and Wiedemann (LRW) [44, 45] 

q L nw = ^^VXT\ (4.2) 

LRW use the dipole formula as a definition of q and compute a strictly lightlike Wilson 
loop. This Wilson loop corresponds into a space-like surface. In this computation special 
care has to be taken in approaching the limits M — > oo and v — > 1. Regardless of the concern 
expressed by some authors about this computation [16, 27, 55] (most of which have been 
addressed in [45]), we do not completely understand the source of this discrepancy. The 



9 There is a trivial factor of 2 difference stemming from the number of spatial dimensions, kt — * \ ■ Further 
q is often expressed in the adjoint representation, so that qA = ^q_F- We will restrict our analysis to the 
fundamental representation and the value of q quoted below differs from the adjoint q in Eq. (15) of 
Ref. [44] by an appropriate factor of two. 
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difference between the two results may lie in the fact that the Wilson line in Ref. [44, 45] is 
strictly light like. 

Certainly the calculation described here is limited to the regime 10 



To see this we recall that the quark moving with velocity v was constructed by slowly turning 
on an electric field and accelerating the quark to its terminal velocity [16]. The equation of 
motion of the quark is 



computed from the Born-Infeld action for the probe brane in the AdS geometry 11 , 



The critical value for the electric field is then 



Equating E ~ rjp we obtain the bound written above, Eq. (4.3). Nevertheless by taking the 
mass to infinity we may approach the lightlike Wilson line from below. 

In summary, the discrepancy may have its origin in the derivation of the dipole formula 
or in the relation between radiative energy loss and the squared momentum transfer. 
Both of these analyses are derived in perturbation theory. A careful analysis of the 
approximations needed to derive these results may resolve the discrepancy and lead to a 
deeper understanding of radiation in a strongly coupled field theory. 

Note Added. During the completion of this work a preprint by S. Gubser, hep-th/0612143, 
appeared and revealed the properties of the world sheet black hole. We gratefully acknowl- 
edge illuminating discussions over the past weeks with Professor Gubser, who pointed out 
(amongst other things) an algebraic error in our draft which had led to a kt independent of 
the quark velocity. Even after correcting this error this manuscript did not agree with the 
first version of hep-th/0612143. The differences have since been resolved. 
Acknowledgments. J.C thanks J. L. Albacete, N. Armesto, V. Koch, C. Salgado, X. N. 
Wang, U. A. Wiedemann for useful discussions about the computation of radiative losses 
and H. Liu about the computation of the light-like Wilson Loops. This work was supported 
by the Director, Office of Science, Office of High Energy and Nuclear Physics, Division of 
Nuclear Physics, and by the Office of Basic Energy Sciences, Division of Nuclear Sciences, 



This paragraph was triggered by discussions with H. Liu, K. Rajagopal and U. Wiedemann. 
We thank H. Liu for pointing out to us the correct value of the critical electric field. 
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